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Abstract 

We study theoretically the topological quantum phase transition in Cavity QED lattice. We 
predict the condition for non-topological phase to the topological phase transition conditions for 
three different model Hamiltonians in cavity QED lattice. We study these topological quantum 
phase transition through winding number, which is a topological invariant quantity. We argue 
that the appearance of topological phase in these systems where the discrete Z 2 symmetry broken. 
We show that the non-topological state is the vacuum state of the system where each cavity 
contains fermionic type excitations from light-matter interaction whereas the topological state of 
system contains Majorana modes of excitations at the end cavity of the lattice. 
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I. INTRODUCTION 


In condensed matter system, quantum phase transition plays a significant role to study 
and explain different quantum phases of the system which explain the different broken 
symmetry states of the system, which are used to describe with the concept of local order 
parameter of Ginzburg-Landau Theory DEI. 

The major limitation of Ginzburg-Landau theory is that the order parameter is treated as 
a local order parameter, this is over come by Topological quantum phase transition [Tj. In 
recent years, it reveals experimentally and theoretically in condensed matter system that 
the typical order parameter is not a local order parameter rather it is highly non-local order 
parameter [1]. One of the example for the non local order parameter is quantum Hall state 
which corresponds to the annihilating an electron at a position by unwinding the number 
of fluxes. This flux unwinding process is highly non-local and Ginzburg-Landau theory is 
not able to describe this phenomena. The states with non-local order parameter are termed 
topologically ordered. The fact that topology is a characterization of global shape and it 
is invariant under a small local deformation it becomes one of the important property of 
the system [1]. This topological robustness protects the quantum state from the external 
perturbation. As a consequence of it several interesting physical phenomena appears in the 
low dimensional quantum many body system Di- 

As the physics of spontaneous symmetry breaking is absent for the topological state of 
matter, the concept of order parameter is also absent to explain the relevant low energy 
physics of these systems. The topological phases are characterized by topological invariants 
integer number. In the present research problem, we do the explicit study of the topological 
quantum phase transition through winding number study. 

Motivation: 

The study of topological phases of quantum many body systems are still in the beginning 
phase. Therefore to search the topological states in different physical system is increasing 
rapidly. The topological phases of condensed matter system attract much more attention 
due to their practical application in low dimensional quantum many body system [1,3]. 
Here we state few examples of topological states of matter in condensed matter physics. 
Integer quantum Hall effect, fractional quantum Hall effect, quantum spin Hall effect, 
topological insulator. But the topological state of matter for cavity QED has not explored 
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yet to the mark HE]. 

The other part of the motivation comes from the Kitaev’s seminal paper [B]. Kitaev has 
proposed the model of one dimensional spinless p-wave superconductors to realize the 
existence of topological phase. At that period, it is difficult to realize Kitaev’s model in 
reality. The electron carry spin-1/2, the first step is to freeze the spin of the particle so that 
the system appears as a one dimensional spinless system. In the interacting light-matter 
system, specially in the cavity QED system where the experimental advancement is in the 
state of art and the spin of the system mimics the state level difference and the measurement 
of quantum state is extremely precise HE]. In low-dimensional interacting light-matter 
system excitations appear as a collective mode under certain physical conditions it behaves 
as a Majorana fermion mode j9]. Therefore we decide to explain the topological state of 
interacting light-matter system and also to realize the Kitaev’s model for cavity QED system. 


II. THE MODEL HAMILTONIANS AND THE DERIVATION OF EFFECTIVE 
HAMILTONIANS 

The Hamiltonian of our present study consists of three parts: 

H = Ha + He + Hac ( 1 ) 

The Hamiltonians are the following 

N 

Ha ^ I T ^ab\^j I: (^) 

i=i 

where j is the cavity index. Uab and Ue are the energies of the state \h > and the excited 
state respectively. The energy level of state \a > is set as zero. \a > and \b > are the two 
stable state of a atom in the cavity and |e > is the excited state of that atom in the same 
cavity. The following Hamiltonian describes the photons in the cavity, 

N N 

He = (Uc^a/aj -h Jc ^(a/aj+i-F h.c), (3) 

i=i j=i 

where ajl(aj) is the photon creation (annihilation) operator for the photon field in the 
j’th cavity, ue is the energy of photons and Je is the tunneling rate of photons between 
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neighboring cavities. The interaction between the atoms and photons and also by the driving 
lasers are described by 

^ Q 

Hac = ga(^j)\Gj >< djl + h.c] + [a -f-)- b]. (4) 

j=i 

Here Qa and gb are the couplings of the cavity mode for the transition from the energy states 
I a > and |6 > to the excited state, fla and fib are the Rabi frequencies of the lasers with 
frequencies uja and Ub respectively. 

The authors of Ref. have derived an effective spin model by considering the following 
physical processes: A virtual process regarding emission and absorption of photons between 
the two stable states of neighboring cavity yields the resulting effective Hamiltonian as 

N N 

Hxy = ^Baj^ + + -^(Tfaj+i- + h.c). ( 5 ) 

i=i j=i 

When J 2 is real then this Hamiltonian reduces to the XY model. Where aj^ = \hj >< 

H - >< = \bj >< = 1 % >< bj\ ■ 

N 

Hxy = + CxAfi+l^) 

i=l 

N 

= ^ ^ + Jy(Ti^CTj+i^). (6) 

i=l 

With Jx = (Ji + J 2 ) and Jy = (Ji — J 2 ). 

We follow the references [7j, to present the analytical expression for the different physical 
parameters of the system. 


Ji = 


72 

4 


\fla?gb^ 




+ 




?ga^ 


Ah 


)) -h 


72 /flaflbgagb 

AaAfe 


(7) 

( 8 ) 


The detail analytical expression for 7 ^ 6 , 71 , 72 , A^, A;,, , 5b^ and fl^ are relegated to the 

appendix. 

Here we discuss very briefly about an effective ^-component of interactions (aAcJi+i^) in such 
a system. The authors of Ref. [Si have proposed the same atomic level conhguration but 
having only one laser of frequency u that mediates the atom-atom coupling through virtual 
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photons. Another laser held with frequency v is used to tune the effective magnetic held. 
In this case the Hamiltonian Hac changes but the Hamiltonians Ha and He are the same. 


O A 

Hac = El(2+ 2e-'''‘‘9aa,)\ej X a,\ + h.c] 
t = l 

[ci 6]. (9) 


Here, and Hb are the Rabi frequencies of the driving laser with frequency oj on transition 
|a>—)■ |e> , |6>—)■ |e>, whereas Ka and Kb are the driving laser with frequency v on 
transition |a>—|e> , |6>—>■ |e>. One can eliminate adiabatically the excited atomic 
levels and photons by considering the interaction picture with respect to Hq = Ha + He 
[6,7]. They have considered the detuning parameter in such a way that the Raman transitions 
between two level are suppressed and also chosen the parameter in such a way that the 
dominant two-photon processes are thus no transition between levels a and b. Whenever 
two atoms exchange a virtual photon both of them experience a Stark shift and play the 
role of an effective interaction [Hi- Then the effective Hamiltonian reduces to 

N 

Hzz = -|- JzCj^aj+i^). ( 10 ) 

i=i 

These two parameters can be tuned independently by varying the laser frequencies. Finally, 
they have obtained an effective model by combining Hamiltonians H^y and H^z by using 
Suzuki-Trotter formalism m- The effective Hamiltonian simulated by this procedure is 

N 


Hspin 

i=i 




( 11 ) 


a=x^y^z 


where Btot = B + B^. It has been shown in Ref. [H] that Jy is less than J^. From the 
analytical expressions of Jx and Jy, it is clear that the magnitudes of Ji and J 2 are different. 

«« = ^ - A) - Al. 

The quantum state engineering of cavity QED is in the state of art due to the rapid progress 
of technological development of this held [3,4]. Therefore one can achieve this limit to get 
the desire Hamiltonian and quantum state of the system, when we consider the situation 
where Jj, = 0 and J 2 = 0. In this limit the atom-photon coupling strength Qa = gt- The 
detail derivation is relegated to the appendix. 
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From the above equation, We get the following relations, gb = ga to get the transverse Ising 
model. The detail derivation is relegated to the appendix. 

N 

z{j) + + ^))- ( 12 ) 

i=i 

One can write the above Hamiltonian in spinless fermion operators by using the Jordan- 
Wigner transformation. To do so, we use the following relation. 

^n+1 (t/’n t/’n)(t/’n+l T t/’n+l) 

an^an+1^ = {'4’j - - ^n+l) 

an^an+l^ = - l)(2^n+lVn+l - 1) (13) 

The Hamiltonian, Hi become. 

Hi = Jx ^(t/’^(n)'^(n + 1) + h.c) + Jx ^(t/’^(n)'^^(n + 1) + h.c) + 2B'^^il)\n)il){n). (14) 

n n n 

We get this Hamiltonian for the condition ga = gb- 

Similarly for the Hamiltonian, iJ 2 , where Jx and Jy are non-zero. One can write the 
Hamiltonian in the following form. 

H 2 = {Jx~\~Jy ) '^{i)\n)i){n+l)+h.c) + \Jx-Jy\y^Xi)\n)i)\n+l)+h.c)+2B E 

n n n 

(15) 

We get this Hamiltonian for the condition, Qb*Qb^a = ^a*ga^b- The detail derivation is 
relegated to the appendix. 

Similarly for the Hamiltonian, H^, where Jx, Jy and are non-zero. 

= {JxJJy) ''^{i)\n)%l){n+l)+h.c) + \Jx-Jy \ ''^{gl)\n)'4}\n+l)+h.c) + {2B+Sp-4:) 

n n n 

(16) 

where p =< > is the density of excitation in interacting light-matter system. Here 

we do the many body physics decoupling scheme to reduce the quartic interaction of the 
Hamiltonian to quadratic one. 

(Jn(yn+i = (8p - A)'4)\n)'4){n) (17) 
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After the Fourier transformation, the Hamiltonian, Hi , reduce to. 


Hi 2 ^ ^ ^2H JxCosk^iyCj^Cj^ A 
k>0 

+2iJx ^ sink{ck^c-k^ + CkC-k). (18) 

A:>0 

Similarly for the Hamiltonian, H 2 reduced to. 


H2 = 25^(25+ (4 + 4 )cosk){ck^Ck + c_fc^c_fc) 

fc >0 

+2^ \Jx '41 ^ ^ StTlk i^Ck C—k d“ CkC—k^ • (^^) 

k>0 

Similarly the Hamiltonian, H^ reduced to. 


H, = 2^((2B + 8p - 4) + (4 + Jy)cosk){ck^Ck + c-k^c-k) 

k >0 

+2i\Jx - 41 y^gw/c(cfc 4 _fc^ + CfcC_fc). (20) 

A:>0 

Now our main interest is to study the topological quantum phase transition in cavity QED 
lattice system based on these models. Our starting point is to recast our three Hamiltonians 
{Hi, H 2 , H 3 ) of interacting light-matter system. 

Hj=ix 3 = -4(c„4n+i + h.c) - - |A|'^(c„c„+i + h.c)]. (21) 

n 

In the above Hamiltonian, we neglect the common negative sign which will not alter the 
relevant physics of the system. 

= 4, = 4 and = 2B. 

t(2) = (j^ + 4 ), a( 2) = (4 - 4) and +(2) = 2B. 
t(3) = (4 + 4), A® = (4 - 4) and +(2) = 2H + 8p - 4. 

The bulk properties of Hamiltonian can be studied in the momentum space. One can write 
down the Hamiltonian in momentum space as. 

-^J=1.2.3 = H^=l,2,3^k (22) 

k 
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where, 


^ / ,(J=1.2,3)(^) 2A(''=l’2.3)*(fc) 

\^2A(^=l’2.3)(;t) -e('^=l’2:3) 

and = —iA^-^^^’^’^^sink. 

This Hamiltonian corresponds to the p-wave snpercondncting phase, one can nnderstand 
this in the following way. 

One can also write down the above Hamiltonian in Bogolnibov energy spectrum, 

^j=i,2,3 = 5^^WVafc (23) 

k 

Here E). is the energy spectrum in bulk and and are the Bogoliubov quasiparticles 
operators. We express the model Hamiltonians of our system in terms of spinless p-wave 
superconducting Hamiltonian thus the starting point of our analysis is the seminal paper of 
Kitaev [ 6 ]. 

In the Majorana fermion basis Hamiltonian reduce to: 

j^J-1,2,3 ^ j ^(1 -)- i'-^B,n'lA,n) + ( 2 ) ^'^)7B,n7A,n+l + (A'^ — )'yA,n'lB,n+l 

Here we use these analytical relations to derive the above Hamiltonian: 

Cji '^B,ii'^A,i} 'lA,i- \Aoi,xi'^l3,y\ ‘^^a,fS^x,y 

The index A and B are the arbitrary index. In the Kitaev’s chain if a Majorana fermion of 
7 a occurs at one end of the chain then the Majorana fermion 7 ^ must occurs at the other 
end of the chain [ 6 ]. 

Now we discuss the non-topological states of the system. The hrst corresponds to /i < 0 but 
t = A = 0. From this condition, we analyze non-topological states of the three Hamiltonians. 
For the Hamiltonian, Hi. 

(1). |f2a| = 0 and \Qb\ = 0, (2). gb = 0,ga = 0, (3). \Qb\ =0,ga = 0. (4). Ifial =0,gb = 0. 

For the Hamiltonian, H 2 , we obtain the same conditions to achieve the non-topological phase 
as we obtain in Hi. But the condition for /i is different for H^. For Hi and H 2 , 61 —2/3 < 0, 
but for H 3 the condition is — 2/3 -|- 8 p — 4 < 0. 

Physical explanation of the non-topological phase is the following: 

The hrst term of the above Hamiltonian yields a coupling between the Majorana fermion 
modes in the same site. In the cavity QED system, this situation corresponds to the different 
kind of light-matter interactions which obeys that the coupled Majorana fermion modes 




t 



t 



FIG. 1: (Color online) Figures show the variation of €b with t for different values A(= d = 1,2, 3,4). 
The upper panel is for the /r = 0, the middle one is for /r = —1.5 and the lower one is for ^ = —2. 
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condition in each cavity. It is well known that the two Majorana fermion modes produce 
a fermion mode. Thus the case of cavity QED system, the non-topological state of the 
system is the fermionic excitations. We consider this state as the vacuum state of the 
system where there is no gapless Majorana fermion excitation states. We consider this non- 
topological state as a vacuum state of the present system, actually it is the conventional 
superconducting phase. 

Now our main intention is to hud out the topological excitation which appears as a Majorana 
fermions. The condition for this phase is p = 0 and t = A ^ 0. Hence the Hamiltonian 
reduced to 

N-l 

H it ^ ^ Ti?,nhA,n-|-l • 

n=l 

It is very clear from the above Hamiltonian that 71 = 7 ^ 1,1 and 72 = 7 _B,Ar are not appear in 
the Hamiltonian. One can also write down the above Hamiltonian by introducing the new 
operator B. 


— 2 (7A,n+l + ilB,n)- 

Af-1 

n=l 

The ends of the chain has zero energy Majorana fermion modes 71 = 'yA,i and 72 = 7b, w- 
These can be considered as a non-local fermion / = |( 7 i + * 72 )- This fermion mode is in 
the zero energy conhguration. In this topological phase, system is in the doubly degenerate 
ground state. One can understand this by following analysis. 

If we consider |0) is a ground state then / |0) = 0 and |1) = |0) is also a ground state 

with opposite fermion parity. The main difference between the conventional superconductor 
and topological superconductor is that the system has a unique ground state with an even 
parity such a way that all the electrons can form cooper pairs. Thus the conventional 
superconducting pairing is the vacuum state of the system. Therefore the non-topological 
state of the system is in the even parity state and the doubly degenerate ground state with 
opposite parity of the system is the topological state of the system. 

Here we consider the most general situation, when p 7 ^ 0 and t ^ A. Here the main intuition 
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is to study the transition for the topological state of the system to the non-topological state 
of the system. 

In this arbitrary limit, the Majorana zero modes are no longer In this limit wave 

function decay exponentially into the bulk of the system. The decay of this wave function 
results in the splitting of the degeneracy between the states (/ |0) = 0, |1) = /t |0)). One 
can write down the effective Hamiltonian as 

Heff = (24) 

_ 

when Cb (X e «o, L is the length of the system, b' and h" are the Majorana fermion at the left 
end and the right end of the chain respectively, is the smallest of |/n|x+|| and |/n|x_||. 
When is 0 then the coherence length ,^o is oo. At this point, the system transit from 
topological state to the non-topological state of the matter. 

—/i ± \/ijfi — -|- 4| Ap 

^ 2(t+ |A|) 

= Mm[|/n|x+||, |/n|xo||] (25) 

e (X e ~when ^ 0 then oc 1. 

In £g. 1, we present the appearance of non-topological trivial state for different values of 
A(A = 1, 2, 3,4). It reveals from our study as we increase the value of A, the non-topological 
phase occurs for the higher value of t. In the present situation the Majorana fermion modes 
decay exponentially into the bulk of the chain. The overlap of these wave function result 

_ L 

in the splitting in the degeneracy between the state |0) and |1) by energy scale e «. This 
figures panel consist of three figures, it is clear from these study that as we go away from 
zero chemical potential to the higher one the peak at t = A of e?, study gradually decreases. 
The coherence length ^ ^ oo for only /x = 0, where the system shows the topological to 
non-topological transition. It is also clear from the analytical expression for et that as we 
increase the length of the system keeping the other physical parameters of the system fixed, 
the system is in the stable topological state. Therefore it is clear from our study that the 
Majorana fermion modes appear at the edge is more stable for the array of larger length 
scale compare to the shorter one. Now we explain the corresponding physics in the light of 
cavity-QED lattice. In the topological phase, there is no bound between Majorana fermion 
mode between the two ends of the cavity. In the case of non-topological state, Majorana 


11 



fermion mode excitation are now bound at the two end of the cavity. It is also clear from the 
above analysis of Kitaev’s formula [H] that the that the transition from topological state to 
non-topological occurs only for /i = 0. (t 7 ^ 0 7 ^ A). But we will also study the topological 
quantum phase transition through the variation of winding number calculation in the next 
section which yields more new and important result. 


TOPOLOGICAL PHASE TRANSITION FROM THE ANALYSIS OF WINDING 
NUMBER 

In this section, we explicitly discuss the physics of topological to non-topological tran¬ 
sition from the analysis of winding number calculation. We already found the analytical 
expressions of three Hamiltonians in momentum space. These Hamiltonians are alike to 
BdG Hamiltonian. 






(26) 


where cos{k), (k) = —2A^'^^sin(k). 

The analytical expression for and A'^(/c) are given in the previous section. 

Topological phase transition can be ascribed by the topological invariant quality. It is 
convenient to dehne this invariant quantity using the Anderson pseudo-spin approach na. 


X(k) ' = l^{kYy + (ex' - 


(27) 


One can write the Hamiltonian as .r where f are Pauli matrices which 

act in the particle-hole basis. It is very clear from the analytical expression that the pseudo 
spin dehne in the YZ plane. 


X{k) 


X{k) 

\x(k)\ 


cos{6k)y + sin{0k)z 


(28) 


Here the momentum states with periodic boundary condition for a ring and the unit 
value x{k) exists on a unit circle in the YZ plane. Therefore 9{k) is a mapping. 
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FIG. 2: ( Color Online), Figures show the variation of winding number with t for different values 
of li. The upper panel is for the /r = 0, the middle panel is for // = 0.5 and the lower panel is for 
/i = 1. Here A = 1. 

and the topological invariant is simply the fnndamental gronp of the maping which is 
jnst the integer winding nnmber. The physical interpretation of this qnantity is that the 
nnit vector x(fc) relates in the YZ-plane aronnd the Brillonin zone. It is only an integer 
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FIG. 3: ( Color Online), Figures show the variation of winding number with t for different values 
of li. The upper panel is for the /r = 0, the middle panel is for // = 0.5 and the lower panel is for 
/r = 1. Here A = 10“^. 
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number and therefore can not vary with smooth deformation of the Hamiltonian so large on 
the quasiparticle gap remains hnite. At the point of topological phase transition the winding 
number changes discontinuously. A topological invariant for x(fc) is then expressed by [IB] 


W 


dk 

dvr 


1 dxfi 

dk ■ 


(29) 


Here a and /3 are Y and Z two components and is anti-symmetric tensor. 

It is very clear from our study, fig.2, that at /x = 0 and t = 0 system shows only the 
non-topological state of cavity QED lattice. Apart from f = 0, we observe that the winding 
number is always one. As we go further away from the zero chemical potential, /x = 0.5,1, 
it reveals from our study that the transition occurs from the non-topological state {W = 0) 
to topological {W = 1) state. The non-topological state persists for a range of t and 
then follows a transition to the topological state. To the best of our knowledge, this 
explicit study of the winding number calculation is absent in the literature of interacting 
light-matter system. 

It is clear from this study that as we go away from /x = 0 the non-topological state 
persist for a wider range of t. It is related with the following relations for three different 
Hamiltonians. 

For the Hamiltonian H 2 , one can write the condition for the persistent of non-topological 
state is the following. 


For the Hamiltonian H^, 
state is the following. 


r r) n '72/^a9b ^b9a\ 


(30) 


one can write the condition for the persistent of non-topological 


_ 2)J + 8p - 4 = |(^ + 5^)'. (31) 

In hg.3, we study the behavior of winding number with t for smaller value of A = 10“^. 
It is clear to us from this study that the system is in the topological state for the zero 
chemical potential. As we away from the zero chemical potential the system is in the 
non-topological state, i.e. there is no sharp topological phase transition in the system. 
Fluctuation of the winding number is extremely large that one cannot predict about the 
dehnite topological phase transition. It can be understand in the following way, as A —)■ 0, 
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the system is simple a fermionic chain with out any spinless p-wave superconductivity in 
the Hamiltonian. Therefore the bulk gap of the system is absent which implies that system 
has no topological state with two Majorana fermion mode at the edge. 
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APPENDIX 


/S 


4A: 

+7i^vr-(o b)] 


4(A, - Afe) 


- Ib9b^ - ll9a^ 


(32) 


_ 1 1 _ 1 1 1 _ 1 

N LDa.fo-Wfc N ^k (LOa+U}b)/2-U)k N ^k {LJa-\-OJb)/2-U)k 

(5l ^ab (^a ^b)/2, UJe ^a* i^^ab *^l)* ^kf 

Sb OJe ^k (s^ab ^1) 5 ^k H“ Jc COSh. 


^2 = Ej=a,b 




^16 Aft 


4Aft 


jUftl^ jUft 


4(Aa - Aft) Aft 4{Aj - Aft) 


E 


I'^jl A 2 \ I 2 9b 
-■^l3b9j )+'lbb ^ 


(a b)] (33) 


Here 71 ^ , 72 , 'Jaa Ibb JV Xlfc W-CJ*. • 


Here = 0, this condition, yields the following analytical relation between the different 
physical parameters of the system. 


9b^a _ 

9aAli lift 

And also from the condition of Jy = 0, i.e., Ji = J 2 . 


(34) 


9aAa _ 

gfftAft flft 

From the above equation. We get the following relations, gb 
model. 


(35) 


Qa to get the transverse Ising 
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Derivation of Eq.26: 

Following Ref. [HI US] , one can write the order parameter as the sum of singlet and 

triplet {P{k)). The singlet and triplet component satisfy the following relation 'ip{k) = ijj{—k) 
and P{k) = —P{—k). One can write the general expression for order parameter as 


A(/c) = i'ilj{k)ay + iP{k)aya 


(36) 


Therefore for the p-wave superconductor, one can write the mean held Hamiltonian as 


H = / dklY^ ejk)ca(k)^cjk) + + h.c] 


(37) 


The hrst part of the Hamiltonian is the single particle energy and the second part is the 
pairing energy. Here we consider, P{k) = (0, —A(A;),0). Finally the Hamiltonian become 

H = j dkl'^ ea{k)ca{kyca{k) - ^ '^[2iA{k)cJ{k)cJ{-k) - 2iA{k)c„{k)c„{-k) (38) 

a a. 

In the above expression, we use A{k) = —A{—k). Finally we can write this hamiltonian in 
Nambu representation 

H=^j dk^'ijjjHk‘ipo.{k) (39) 

a. 

'ipa{k) = {ca{k),—Ca^). The half factor is to balance the double counting. Hk = 

e„(/c) -2iA^{k) 

2iAa{k) —ea{k) 
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